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Abstract. We present several arguments which favor the scenario of two coexisting families of compact 
stars: hadronic stars and quark stars. Besides the well known hyperon puzzle of the physics of compact 
stars, a similar puzzle exists also when considering delta resonances. We show that these particles appear 
at densities close to twice saturation density and must be therefore included in the calculations of the 
hadronic equation of state. Such an early appearance is strictly related to the value of the L parameter 
of the symmetry energy that has been found, in recent phenomenological studies, to lie in the range 
40 < L < 62 MeV. We discuss also the threshold for the formation of deltas and hyperons for hot and 
lepton rich hadronic matter. Similarly to the case of hyperons, also delta resonances cause a softening of 
the equation of state which makes it difficult to obtain massive hadronic stars. Quark stars, on the other 
hand, can reach masses up to 2.75M© as predicted by perturbative QCD calculations. We then discuss the 
observational constraints on the masses and the radii of compact stars. The tension between the precise 
measurements of high masses and the indications of the existence of very compact stellar objects (with radii 
of the order of 10 km) is relieved when assuming that very massive compact stars are quark stars and very 
compact stars are hadronic stars. Finally, we discuss recent interesting measurements of the eccentricities 
of the orbits of millisecond pulsars in low mass X-ray binaries. The high values of the eccentricities found 
in some cases could be explained by assuming that the hadronic star, initially present in the binary system, 
converts to a quark star due to the increase of its central density. 

PACS. PACS-key describing text of that key - PACS-key di-scribing text of that key 


1 Introduction 

Ultra-relativistic heavy ions experiments have provided 
many indications of the formation of a new phase of strongly 
interacting matter, named quark gluon plasma, which is 
obtained by heating up hadronic matter (with almost van¬ 
ishing baryon density) to temperatures of few hundreds 
MeV [1]. In this state, the fundamental degrees of freedom 
of QCD, quarks and gluons, are deconfined. A very inter¬ 
esting question of nuclear and hadronic physics concerns 
the possibility of the formation of a deconfined phase also 
at large baryon densities and small temperatures. Natu¬ 
ral systems to look for this state of matter are neutron 
stars. In this respect, the recent discoveries of two stellar 
objects m with masses of M = 2Mq are very promis¬ 
ing: the larger the mass the larger the baryon density in 
the core of the stars. Those massive compact stars are 
therefore the best systems for studying the structure of 
the QCD phase diagram at high densities. 

There is huge collection of theoretical and phenomeno¬ 
logical calculations aiming at establishing whether quark 

Send offprint requests to-. 


matter can form in these stellar objects and which would 
be the possible associated observational signatures. One 
can distinguish three possible scenarios: i) there is only 
one family of compact stars which are hybrid stars i.e. 
stars composed by hadronic matter at low density and 
quark matter at high density. See for instance Refs. [U 
[5116117] for recent calculations providing equations of state 
stiff enough to support a star of 2Mq. ii) high mass “twin 
compact stars” [5] . In this scenario most of the stars would 
be composed only of nucleonic matter. Stellar configura¬ 
tions with masses of about 2Mq could be composed either 
of nucleonic or of hybrid matter. The two coexisting stellar 
conhgurations have different radii: hybrid stars are more 
compact than their “twin hadronic stars”. iii) two sepa¬ 
rated families of compact stars [ 3 : hadronic stars which 
can be very compact and have a maximum mass of about 
1.5 — I. 6 M 0 ; quark stars which can be very massive, up 
to 2 . 75 M 0 , and have large radii. 

The third scenario is the one we will present in this 
paper and in the accompanying paper 2. We will discuss 
in particular the phenomenological motivations in favor of 
this model and we will try to analyze the predictions which 
distinguish this model from the other two. The paper is 
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organized as follows: in Sec.2 we present the calculations 
of the equations of state of hadronic and quark matter. In 
Sec.3 we compare our theoretical results with the observa¬ 
tional constraints on the masses and the radii of compact 
stars. In Sec. 4 we discuss the information one can obtain 
on the Equation of State (EoS) from the study of compact 
stars in binary systems. Einal discussions and conclusions 
are presented in Sec. 5. 


2 Equations of state 

In this Section we present the EoS for hadronic and quark 
matter used in the present investigation. Eor hadronic 
matter a problem hugely discussed in the literature con¬ 
cerns how to reconcile the unavoidable appearance of hy- 
perons (at densities of the order of 2 — 3no) with the exis¬ 
tence of compact stars with masses of 2Mq. This problem 
is discussed in the contributions of Chatterjee and Vidana 
and of Oertel et al. in this volume. The similar problem 
of the appearance of delta resonances at finite density is 
much less discussed in the literature for reasons that we 
will clarify in the following. Here we summarize our find¬ 
ings concerning the simultaneous formation of hyperons 
and deltas within a relativistic mean field approach. 

In the first subsection, the /3-stable hadronic EoS is 
studied in the regime of zero temperature within a nonlin¬ 
ear relativistic Walecka-type model and in the framework 
of the so-called SEHo parametrization which takes into 
account the recent experimental constraints. In the sec¬ 
ond subsection, we extend the study of the SEHo model 
at finite value of entropy per baryon. 


density dependence of the symmetry energy at supranu¬ 
clear densities and satisfy all of the experimental con¬ 
straints both from terrestrial and astrophysical data by 
restricting L to the range of 40 MeV ^ L < 62 MeV [TT] 
Q. To this purpose, we are going to compare the results 
in the framework of the GM3 model (with a value of L ~ 
80 MeV, automatically fixed once a specific value of S is 
adopted) with a more sophisticated EoS, called SEHo, for 
which 5 = 32 MeV (very close to the GM3 value) and L 
= 47 MeV 

In the GM3 model the general form of lagrangian is 
given by m 

jC-octet = '^^k -mk+ gcrkO- - 

k 

gpkit^Y' ~ 

-k 

+ + U{a,u},p), (2) 

where the index k runs over the baryon octet, mk is the 
bare mass of the baryon k, Tk is the isospin operator and 
finally U is the mesons potential which can contain non 
linear interaction terms. 

Goncerning hyperons, with the exception of the A, 
their binding energies in hypernuclei are highly uncer¬ 
tain (see, for example. Ref. [TS] and references therein) 
and thus also their couplings with mesons are poorly con¬ 
strained. Here, we use the parameters set of Refs. [T6l 
[nii obtained by reproducing the following values of the 
binding energies in nuclear matter Uf 


2.1 Hadronic equations of state including Z\-isobars 
and hyperons at T = 0 


= -28 MeV, C/f = 30 MeV, US = -18 MeV. (3) 

For the coupling with vector mesons we use the SU(6) 
symmetry relations: 


Concerning the hadronic EoS we consider in this paper 
two different relativistic EoS with the inclusion of the 
octet of lightest baryons (nucleons and hyperons) and Z\(1232) 
isobar resonances. First, we study the nonlinear GM3 model 
of Glendenning-Moszkowsky in which the interaction be¬ 
tween baryons is mediated by the exchange of a scalar 
meson ct, an isoscalar vector meson w and a isovector vec¬ 
tor p [To]. Let us note that, within the GM3 parametriza¬ 
tion, only the experimental value of the symmetry energy 
at saturation S is used to fix the coupling between the p 
meson and the nucleons. However, recently, a remarkable 
concordance among experimental, theoretical, and obser¬ 
vational studies has been found im, by allowing to signif¬ 
icantly constrain also the value of L, the derivative with 
respect to the density of the symmetry energy S at satu¬ 
ration: 

d*5* 

L = “ins \nB=no ■ ( 1 ) 


1 

1 

1 

(4) 

-^guN = 


— — 9u)3 

9pN = 

^9pE 

— gujE 1 9pA — 0 ■ 

(5) 


In relativistic heavy ion collisions, where large values 
of temperature and density can be reached, a state of reso¬ 
nance matter may be formed and the Z\(I232)-isobars are 
expected to play a central role. [T51[T^[^II2 1 ll^ . More¬ 
over, it has been pointed out that the existence of Z\s can 
be very relevant also in the core of neutron stars cniEii 

miMm- 

The mean-field Lagrangian density for the Z\-isobars 
can be then expressed as 

CA=i’Au - 9crAO-) “ gujAlp.^^^ 

-gpAip.hp'^] i^A > ( 6 ) 


Therefore, extensions of the GM relativistic mean-field 
model have been implemented which include p meson self¬ 
interaction terms. These new parametrization modify the 


^ Notice that this rather restricted range for the values of L is 
still debated, see for instance Ref. m where a looser constraint 
is proposed. 
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where is the Rarita-Schwinger spinor for the Z\-isobars 
{A++, Z\+, A°, A-) and h = diag(3/2,1/2, -1/2, -3/2) 
is the matrix containing the isospin charges of the Z\s. 

As customary, for the couplings of hyperons and A 
isobars with the mesons, we introduce the ratios 

^ai — 9<7 i / 9<yN 5 i j9^N ^pi — 9pi!9pN ^ ( 7 ) 

where the index i runs over all the hyperons and A isobars. 

Concerning the values of the Z\-meson couplings, if the 
SU(6) symmetry is exact, one adopts the universal cou¬ 
plings Xa-A = XcjA = 1. As already extensively discussed 
in Ref. [55], among the four A isobars, the is likely 
to appear first because it can replace a neutron and an 
electron at the top of their Fermi seas in /3-stable matter. 
However, this particle is “isospin unfavored” because its 
isospin charge = —3/2 has the same sign of the isospin 
charge of the neutron. For large values of the symmetry 
energy S and, therefore, of 9pA: the A~ appears at very 
large densities or it does not appear at all in dense matter 
thus playing no role in compact stars. Indeed, in Ref. [55] 
the Z\-isobars could appear in neutron stars only for not 
physical small values of the symmetry energy, obtained 
by setting gpi = 0 for all the baryons. However, as al¬ 
ready observed, in the GM3 model the coupling gppf is 
fixed by using the experimental value of the symmetry 
energy, the most recent estimates ranging in the interval 
29<S< 32.7 MeV m- In this scheme no experimental 
information on the density dependence of the symmetry 
energy can be incorporated and in particular the L pa¬ 
rameter is automatically fixed once a specific value of S 
is adopted. It turns out that in the models introduced 
in Refs. mm. L ~ 80 MeV and it is thus significantly 
higher than the values suggested by the most recent anal¬ 
ysis [TT] . 

Moreover, in this context let us observe that the SU(6) 
symmetry is not exactly fulfilled and one may assume 
the scalar coupling ratio x^a > 1 with a value close to 
the mass ratio of the A and the nucleon [5^. On the 
other hand, QCD finite-density sum rule results show that 
the Lorentz vector self-energy for the A is significantly 
smaller than the nucleon vector self-energy implying there¬ 
fore XuiA < 1 [30] ■ 

In the many body analysis of Ref. [31] , the real part of 
the A self-energy has been evaluated to be about —30 MeV 
at hb = 0.75 riQ. Notice that this self energy is relative to 
the one of the nucleon and the total potential felt by the A 
is the sum of its self energy and of the nucleon potential, 
a number of the order of —80 MeV. Also phenomenolog¬ 
ical analysis have been performed of data from electron- 
nucleus [35 1 I53 1 IM]. photo-absorption [5^ and pion-nucleus 
scattering [551157] . Such analyses suggest a more attractive 
interaction of the A in the nuclear medium with respect to 
the nucleon one (see Refs. mm for more details). New 
analyses, and possibly new experiments, aiming at a bet¬ 
ter determinations of these couplings would be extremely 
important. Notice also that no information is available for 
XpA which in principle could be extracted by analyzing 
scattering on neutron rich nuclei (see the recent discus¬ 
sion in gnus])- 



Fig. 1. Particles fractions as functions of baryon density (in 
units of the nuclear saturation density no) in the GM3 model 
for Xa-A = 1.25, XuiA — 1. 


The threshold for the formation of the i-th baryon is 
given by the following relation: 

9-t>mi- g„ia + gu,^uJ + higpiP , ( 8 ) 

where cr, w and p are the expectation values of the cor¬ 
responding fields, /ii, rrii and t^i are the chemical po¬ 
tential, the mass and the isospin charge of the baryons. 
The baryon chemical potential pi are obtained by the /3- 
equilibrium conditions: 

Pi = PB+ Ci PC , (9) 

where pB and pc are the chemical potentials associated 
with the conservation of the baryon number and the elec¬ 
tric charge respectively and Ci is the electric charge of the 
i-th baryon. 

In Fig. 1, we display the baryon density dependence 
of the particle’s fractions in the GM3 model for XrrA = 
1.25, Xi^A = 1 and neglecting in this case the coupling of 
the Z\-isobars with the p meson (xpA = 0). Let us note 
that in this scheme the appearance of the Z\-isobars is 
a consequence of the introduction of a more attractive 
interaction [x^a > 1) of Z\-particles with respect to the 
nucleon in the mean field approximation, as in Refs. |291 
[5U1I55] . It is remarkable that the early appearance of A 
resonances, the first one being the Z\“, considerably shifts 
the onset of hyperons which start to form at densities of 
^ 5 po (see the curve for the A’s). 

With the purpose of including the new experimental 
constraints about the value of the density derivative of 
the symmetry energy L, we can first study an extended 
GM3 model by considering the following density depen¬ 
dent baryon-p meson coupling [551155] 

= ( 10 ) 

In this way we introduce a single parameter a which affects 
only the value of L leaving untouched the other properties 
of nuclear matter at saturation. 
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The role of such density dependent baryon-p meson 
coupling in the modified GM3 model can be observed in 
Fig. [21 where is reported the value of for the different 
baryons as a function of L. We limit this first discussion 
to the case of the A, and S~ which are the first heavy 
baryons appearing as the density increases (notice that E 
hyperons are unfavored due to their repulsive potential) 
in the so-called universal coupling x^a = XujA = XpA = 1- 
One can notice the different behavior of the thresholds: 
the larger the value of L the larger and the smaller 

niit and 

At high values of L, larger than about 65 MeV, the 
threshold of the A~ increases very rapidly with L. This 
corresponds to the values of L for which the appears 
before the A~ thus completely suppressing those particles. 
Indeed within the GM3 model, for which L ~ 80 MeV, the 
Z\“ do not appear at all as already found in Ref. [55] • Sim¬ 
ilarly, one can notice that if the isobars are formed before 
the hyperons, what happens below L ^ 56 MeV, and 
are shifted to larger densities, as already noticed in 
Ref. (9). Analogues results have been found in Ref. [28] . 
where two cases are analyzed, corresponding to a finite 
and to a vanishing value of gpN, with the result that in 
the case of gpN = 0 the isobars are favored. The blue lines 
mark the range of the values of L indicated by the analysis 
of Ref. m- Therefore, the recent constraints on L imply 
that at densities close to three times no both the hyperons 
and the isobars must be included in the equation of state 
and for the lower allowed values of L, the isobars appear 
even before the hyperons. Finally, let us stress that in this 
analysis we have chosen a rather conservative choice for 
the couplings between Z\s and mesons. If higher values of 
XaA and or lower values for x^iA are adopted, can re¬ 
sult to be smaller than and n“ for all the acceptable 
values of L. 

In comparison with the previous results, we consider 
now a more sophisticated model for the EoS proposed >h' 
in Ref. [T5irn] . where we use the parametrization called 
SFHo for which S' = 32 MeV (very close to the GM3 
value) and L = 47 MeV with the addition of hyperons 
(assuming SU(6) symmetry) and Z\-isobar degrees of free¬ 
dom (assuming Xa-A = XpA = 1 and different values for 

XujA ) ■ 

Results for the particles’ fractions as function of the 
baryon density in /3-stable matter within the SFHo model 
are displayed in Fig. [3] In the upper panel, we have in¬ 
cluded only hyperons: the A and the S~ appear at a den¬ 
sity of about 0.5 fm“^ and then the 5"° at a density of 
about 1.1 fm“^. In the lower panel we include also the 
A isobars. In agreement with the hndings of the previous 
analysis, for small values of L the As appear at densities 
relevant for neutron stars and actually, in the SFHo model, 
they appear even before the hyperons with the A~ formed 
at a density of about 0.4 fm“^. The appearance of these 
particles delays the appearance of hyperons. It is impor¬ 
tant to remark that, within the SFHo model, even using 
more repulsive interaction than nucleons, x^jA = 1-1, the 
A~ appear before hyperons. 



Fig. 2. Threshold densities of hyperons and As as functions 
of the L parameter for XaA = XuiA = XpA ~ 1 within the mod¬ 
ified GM3 model with the introduction of a density dependent 
baryon-p meson coupling (see the text for details). The contin¬ 
uous lines refer to the case in which all the degrees of freedom 
are included in the computation of the EoS and the dashed 
lines refer to the case in which either hyperons or As are ar¬ 
tificially switched off. The vertical lines indicate the range of 
allowed values of L, as found in m 



Fig. 3. Particles fractions as functions of the baryon density 
within the SFHo model: only hyperons (panel (a)), hyperons 
and As (panel (b)) for Xo-a = x^ja = XpA = 1- The red line 
indicates the fraction of the A~ which among the four As 
are the first to appear. The blue and the green vertical lines 
indicate the onset of the formations of A~ for Xc^a = 0.9 and 
XojA = 1.1, respectively. 


In conclusion, the early appearance of A-isobars re¬ 
sults to be strictly related to the value of the L parameter 
of the symmetry energy and, for the values of L indicated 
in Ref. m. such particle degrees of freedom influence the 
appearance of hyperons and cannot be neglected in the 
EoS. These results have been confirmed in the more re¬ 
cent Ref. m- 
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2.2 Hadronic equation of state at finite entropy per 
baryon 


In this subsection we are going to study the behavior of 
the hadronic EoS for conditions realized in protoneutron 
stars (PNS). In particular, we focus our investigation by 
considering the more realistic SFHo parametrization in 
the first stage of the protoneutron stars (PNS) evolution, 
corresponding to a total entropy per baryon equal to one, 
in which neutrinos are trapped and strongly influence the 
chemical composition of the PNS. Therefore, we also take 
into account of leptons particle by fixing the lepton frac¬ 
tion 

Pl = Pe + , (11) 


where rie, and ns are the electron, neutrino and baryon 
number densities, respectively. 

The total entropy per baryon is calculated by means 


of 


Sb+Si 

® - 

T pb 


( 12 ) 


where Sb = Pb + sb - Y.i=B Si = Pi + ei - 

Si=/ PiPii ctnd the sums are extended over all the baryons 
and leptons species. 

It is well known that the presence of trapped neutrinos 
significantly alter the protons and the electrons abundance 
and strongly influence the threshold of hyperons forma¬ 
tion. This is also true in the presence of Zl-isobar degrees 
of freedom. With the purpose of investigating this prob¬ 
lem, in Fig. |4j we report the particle concentrations Yi as 
a function of the baryon density for s = 1 and Yb = 0.4 
in the SFHo parametrization with the coupling XrrA = 1-0 
(upper panel) and XcrA = 1-1 (lower panel). 

Let us observe that, in the case of the universal cou¬ 
pling XaA = Xi^A = XpA = 1 (upper panel), A and A~ 
particles appear approximately at the same baryon den¬ 
sity {ub ~ 3no). For x^a = 1-1 (lower panel), the onset 
of A~ particles is shifted at lower densities and the pres¬ 
ence of Z\-isobars become more relevant. In both cases, the 
population of strange A particle becomes relevant (greater 
than 5%) at about ub ~ 4no- 

The features observed in the particle concentration are 
also reflected in Fig. [SI where we show the temperature 
as a function of the baryon density for nucleonic matter 
(np), for hyperonic matter (npH) and with the inclusion 
of Z\-isobar degrees of freedom (npHA). As before, in the 
presence of A particles, we have considered two different 
meson-Z\ couplings (continuous line for x^a = 1-0 and 
dashed line for XaA = 1-1). For x^a = 1-0, A and A~ 
particles start at ng ss 3 no and, for a large baryon density 
range, the behavior is almost isothermal (T w 18 A 20 
MeV). A more discontinuous behavior can be observed in 
the case of x^a = 1-1, due to the presence of the four 
A-isobar states. 

In Fig. |BJ the gravitational mass as a function of the 
central baryon density ric for s = 1 and Yb = 0.4 is re¬ 
ported for two different values of the Xo-a coupling ra¬ 
tio. In this case no appreciable difference can be observed 
for npH and npHA curves (overlapped blue and green 
curves in the figure), except for a greater central density 


Yi 



Yi 



Fig. 4. Particle concentrations Yi as a function of the baryon 
density for s = 1 and Yl = 0.4 in the SFHo parametrization 
with the coupling x^a = 1-0 (upper panel) and x^a = 1-1 
(lower panel). In both cases Xu;A = XpA ~ 1- 


ric reached in presence of A particles. In agreement with 
the previous results, we can see from the figure that hy¬ 
perons and As appear at ub ~ 3 no, corresponding to 
Mg ~ 1.45Mq. For stellar configurations with masses be¬ 
low this value, deltas and hyperons do not appear or play 
a marginal role. Therefore we do not expect any difference 
concerning the SN explosion mechanism in the two fami¬ 
lies scenario with respect to the standard one. For larger 
masses there are a few possibilities: 

— hyperons appear and trigger the transition to quark 
matter halting the collapse; 

— hyperons appear but their abundance is not large enough 
to trigger the conversion and a black hole will form af¬ 
ter deleptonization (4^. 

In this little discussion we have not taken into account the 
rotation of the progenitor what can play an important role 
as discussed in paper 2. 


2.3 Quark matter equation of state 

The quark matter EoS at densities reachable in the core 
of compact stars is basically completely unknown. In the 
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T [MeV] 



Fig. 5. Temperature as a function of the baryon density (in 
units of no) at fixed entropy per baryon and lepton fraction 
(s = 1, Vl = 0.4) for the SFHo parametrization. The labels np, 
npH and npHA stand for nucleons, nucleons plus hyperons, 
nucleons plus hyperons and A-isobars, respectively. The results 
for two different meson-A couplings in the curves npHA are 
reported, the continuous line refers to x^a = 1.0 while the 
dashed line stands for x^a = 1.1. 


Mg/Mo 



Fig. 6. Gravitational mass as a function of the central baryon 
density ric for s = 1 and Yl = 0.4 in the SFHo parametrization 
with the couplings XaA = 1.0 and x^a ~ 1.1. 


literature, bag models or chiral models at finite chemi¬ 
cal potential have been widely used which capture two 
important non-perturbative aspects of QCD: confinement 
the former and chiral symmetry breaking the latter m- 
Alternatively, one can resort to perturbative QCD cal¬ 
culations: it has been shown that at high temperature, 
vanishing chemical potential and finite quark masses per¬ 
turbative calculations provide results consistent with lat¬ 
tice QCD for temperature larger than about 0.2 GeV SZl 
EH]. The same technique has been used to compute the 
EoS at finite chemical potential and vanishing tempera¬ 
ture with the very interesting result that the strange quark 
matter EoS could be stiff enough to support stars with 
2.75Mq [17]. Here, we adopt the parametrization of the 
pQCD calculations at finite chemical potential presented 


in [48] . This parametrization has only one free parame¬ 
ter, the scale parameter X, which is the ratio between the 
renormalization scale and the baryon chemical potential 
and it ranges between 1 and 4. We choose here the value 
X = 3.5 for which the maximum mass of quark stars 
is 2.53M0 (see solid red line in FigjTj). It is important 
to remark that all the calculations of the quark matter 
EoS, such as the ones of m. are affected by large un¬ 
certainties which inevitably affect also the predicted value 
of the maximum mass. For instance, also in the simple 
MIT bag model parametrization of Ref. [49] . it has been 
proven that very massive configurations can be obtained 
when varying the free parameters of the model. Moreover, 
it has been shown that, color superconductivity, which is a 
purely non-perturbative phenomenon, helps in obtaining 
high masses. Recently, within a SU(3) quark meson model 
it has also been shown that massive quark stars configura¬ 
tions are allowed although in a specific parameters region 
m- Clearly, new and precise astrophysical measurements 
are needed to constrain the properties of dense matter. 

3 Masses and radii: theory vs observations 

Let us discuss the information on the mass-radius rela¬ 
tion obtained by means of astrophysical observations. The 
direct and precise measurements of the masses of PSR 
J1614-2230 with M = 1.97 ± 0.04Mo |2] and of PSR 
J0348-I-0432, with M = 2.01 ± O.O4M0 |5], clearly repre¬ 
sent the most important constraints that theoretical cal¬ 
culations must fulfill. A possible candidate with a mass 
even larger already exists: it is the black widow pulsar 
PSR B1957-I-20 whose mass could be of 2.4 ± O.12M0 
provided that the modeling of the light curves is correct 
m- Taking into account the systematic uncertainties on 
the light curves fit, the lowest limit for its mass turns 
out to be of 1.66M0. Other hints for the existence of 
compact stars heavier than 2Mq have been obtained also 
from the observation and the modeling of short-gamma- 
ray bursts (GRB). The SWIFT experiment has detected 
tens of short-GRB whose light curves display extended 
emissions. X-ray flares and internal plateaus with rapid 
decay at the end of the plateaus (see [52] and Refs, therein). 
These observations favor a model for the inner engine of 
these events which is based on a rapidly spinning mil¬ 
lisecond magnetar formed from the merger of two neutron 
stars. Interestingly, the same stellar objects, but formed 
after a supernova, could be the inner engine of long GRBs 
[55] (see paper 2). In [551155] . the detailed modeling of 
the plateaus seen in short GRBs has provided an impor¬ 
tant constraint on the peak of the expected mass distribu¬ 
tion for the post-merger remnant: M = 2.4:6^q\^Mq (dis¬ 
played in Fig. [7]). Although this limit is not obtained via 
a direct mass measurement and it includes supramassive 
stars, it represents a strong indication that the maximum 
mass of compact stars is significantly larger than 2Mq. 

Let us discuss now radii measurements. One has to re¬ 
mark that radii measurements are much more uncertain 
than mass measurements and all the observational con¬ 
straints are based on specific assumptions made for mod- 
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eling the spectra of the X-ray emissions. In Refs, 
the fits on the thermal emission of 6 quiescent low-mass 
X-ray binaries, under the assumption that all of them have 
the same radius i?, provide the constraint R = 9.4 ± 1.2 
km. We remark however that these results are under de¬ 
bate, see Refs. [S71ISE]. Other indications of the existence 
of stars with small radii can be found from the analysis of 
X-ray bursts in Refs. |591IBD] . In particular, at Icr, the anal¬ 
ysis of [35] indicates radii of about 9.5 ± 1.5km and masses 
of about M = 1 . 6 ±O. 2 M 0 (a previous analysis of 4U1820- 
30 presented in Ref. m has also found rather small radii: 
R = 11.2+°;4km and M = 1.29ti^;J^MQ ). For Cyg X-2 
a radius of about 9 ± 0.5km is inferred for the canonical 
mass of 1.44 ± O.OGMq while for 4U 1728-34 the sug¬ 
gested ranges are 8.7 — 9.7 km for radius and 1.2 — 1.6Mq 
for the mass m- Also, in the analysis of the X-ray pul¬ 
sations of SAX J1808.4-3658 one obtains indications of a 
small radius although the uncertainty is still quite large: 
at 3 ( 7 , O. 8 M 0 < M < 1.7Mq and 5km < R < 13km [531 
153] . These data and the mass-radius relation of hadronic 
stars (with deltas and hyperons) are displayed in FiglS] 

On the other hand, significantly larger radii are ob¬ 
tained by means of pulse phase-resolved X-ray spectroscopy 
of PSR J0437-4715 [53]: the radius is constrained to be 
larger than ~ 14 km at Icr conhdence level assuming the 
mass of the star to be of 1.76M0 (this value has been ob¬ 
tained via the radio timing technique in [55]). Similarly, in 
Ref. [57] , a radius larger than about 14 km is obtained for 
the system RX J1856.53754 by assuming a mass between 
1.5 and I. 8 M 0 . 

The existence of very massive compact stars and the 
possibility that some neutron stars are very compact rep¬ 
resents a serious problem for the theoretical modeling of 
the EoS of strongly interacting matter. While massive 
compact stars imply that the EoS is stiff, a soft EoS is 
instead needed to obtain small radii. This tension is re¬ 
lieved, as proposed in [5], if one assumes that there are 
two families of coexisting compact stars: hadronic stars 
which can be very compact and quark stars which can be 
very massive. Specifically, the massive PSR J1614-2230 
and PSR J0348-I-0432, are interpreted in our scenario as 
quark stars. Similarly, stars with large radii, as the ones 
inferred in the analysis of Refs. [531157] . are again inter¬ 
preted as quark stars. On the other hand, the compact 
stellar objects such as the ones discusses in the analysis of 
[331[351[351l6(]ll611153] would be instead hadronic stars. The 
idea that the equation of state for dense matter has a two- 
phase nature allowing both large and small compact stars 
has been suggested also in the analysis of Ref. [ 68 ] . 

In Fig. [7] and [5] we display the observational constraints 
discussed above and three examples of theoretical mass- 
radius relations (solid and dashed lines) based on the SFHo 
model with A only, or with both A and hyperons for the 
hadronic EoS (here Xa-A = 1.15) and the parametriza- 
tion presented in [33] for the quark matter EoS. All the 
constraints are fulfilled if one assumes that hadronic stars 
and quark stars coexist. Notice that hadronic stars cannot 
reach masses larger than about 1.5 — I. 6 M 0 as a result of 
the softening caused by the formation of deltas and hy¬ 


perons. The so called hyperon puzzle, and similarly the 
delta puzzle pointed out in [38] , is easily solved in our two 
families scenario: hyperons and deltas do reduce the maxi¬ 
mum mass of compact stars to values significantly smaller 
than 2Mq but this fact does not represent a puzzle since 
the most massive objects are actually quark stars. Notice 
that also in the scenarios i) and ii) discussed in the in¬ 
troduction, massive stars are composed mostly of quark 
matter. 

A natural question concerns the way the quark star 
branch is populated. The stellar configuration for which 
the solid black line starts to deviate from the solid green 
line corresponds to the onset of hyperons. Once a critical 
amount of hyperons is present in the center of the star, 
nucleation of quark matter can start and can subsequently 
trigger the conversion to a quark star. The conversion oc¬ 
curs because it is energetically convenient: at a fixed value 
of the baryonic mass, the gravitational mass of the star 
on the quark star branch is smaller than the one on the 
hadronic star branch (see the brown line for an example 
of this conversion process). The process of conversion is 
specifically analyzed in the accompanying paper 2 . 

4 Binary systems 

Another useful constraint on the EoS can be obtained also 
by studying the double pulsar system J0737-3039 [69]. The 
mass of the so called Pulsar B is of 1.249 ± O.OOIM 0 . Un¬ 
der the assumption that this pulsar was formed from an 
electron capture supernova one can infer a baryonic mass 
in the range 1.366 — 1.375Mq [75]. In our scenario we can 
interpret this star as a hadronic star (which contains A 
resonances but which is too light to allow the formation 
of hyperons). In FigiS] we display the relation between 
the gravitational mass and the baryonic mass for hadronic 
stars and quark stars. In the insert, we show also the con¬ 
straint of [75]. Our hadronic EoS is perfectly in agreement 
with the analysis of Ref. [75]: the appearance of delta reso¬ 
nances gives a small additional contribution to the binding 
energy of hadronic stars as compared to nucleonic stars. 
When considering only nucleons (see green line), the con¬ 
straint of [75] is not fulfilled. Notice however that detailed 
supernova simulations have shown that the uncertainties 
associated with the EoS and the wind ablation are such 
that the allowed baryonic mass window is shifted towards 
smaller values m- 

Low-mass X-ray binaries offer another possible hint 
for the existence of two families of compact stars. These 
systems are most probably at the origin of millisecond 
pulsars: within the so called recycling scenario, the neu¬ 
tron star is spun up to milliseconds periods due to the 
accretion of mass from its white dwarf companion. Tidal 
interactions during the accretion phase are responsible for 
the circularization of the orbit and indeed most of the mil¬ 
lisecond pulsars are in circular orbits (with eccentricity e 
from 10“^ to 10“^). However, recently, few examples have 
been discovered having a much larger eccentricity such 
as PSR J2234-I-06 for which e = 0.13 (see other exam¬ 
ples in m)- In this system the white dwarf has a mass 
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of 0.23Mq. The existence of systems with high eccentric¬ 
ities represents a puzzle in the recycling model of pul¬ 
sars. A possible explanation is that the accreting object, 
at some point during its evolution, collapses to a more sta¬ 
ble configuration thus increasing abruptly the eccentricity 
of the binary. In Ref. m it has been investigated the sce¬ 
nario of a rotationally-delayed accretion- induced collapse 
of a super-Chandrasekhar mass white dwarf. In Ref. m 
instead, the accreting star is a neutron star which, due to 
mass accretion, converts into a quark star. In our two fam¬ 
ilies scenario the conversion of a hadronic star to a quark 
star is necessary once a sufficient amount of strangeness is 
formed at the center of the star. As one can see in Figs 
the conversion would occur for masses of the hadronic star 
between 1.35 — 1.6Mq with an energy released in the 
conversion given by the difference between the gravita¬ 
tional mass of the hadronic star Mh and the gravitational 
mass of the quark star Mq computed at the same fixed 
baryonic mass: AM = Mh — Mq ^ O.15M0 (see Fig. IH]). 
The eccentricity is related to the masses of the hadronic 
star, of the quark star and of the white dwarf companion 
Mwd by the following relation: e = AM/{Mq -|- Mwd)- 
It results that e ^ 0.1 if one takes Mwd = 0.23Mq, 
Mh = I. 55 M 0 Mq = 1.4M0. A correction to the eccen¬ 
tricity of the order of ±0.03 is obtained when considering 
that during the conversion the newly born stellar object 
could get a small kick velocity Vk of the order of Ikm/sec 
m- These simple estimates show that in our model the 
values of the eccentricities are quite close to the measured 
ones. It would be therefore interesting to investigate more 
in detail this problem. Future measurements of the masses 
of the compact stars in those eccentric systems will be cru¬ 
cial to test our scenario. 



R[km] 


Fig. 7. Gravitational mass-radius (solid lines) and baryonic 
mass-radius (dashed lines) relations for hadronic stars (includ¬ 
ing only deltas and hyperons and deltas ) and quark stars. 
Some of the most recent observational constraints are also dis¬ 
played (see text). The brown lines show that given a hadronic 
star configuration (in which also hyperons are present), the 
quark star with the same baryon mass has a smaller gravita¬ 
tional mass even if its radius is larger. The conversion of the 
hadronic star into a quark star would be therefore energetically 
favored. 


-- M^-SFHo-HA 


5 Discussion and conclusions 


Ozel et al. 2010 


We have discussed several hints of the existence of two 
coexisting families of compact stars, hadronic stars and 
quark stars. A first important and widely discussed argu¬ 
ment in favor of this scenario is the necessary appearance 
of delta resonances and hyperons as the central density of 
a hadronic star reaches values larger than about 2no. The 
formation of these particles softens the EoS and reduces 
the maximum mass with respect to stars made only of 
nucleons. How much the maximum mass is reduced due 
to the appearance of these particles is the subject of a 
lively and on-going research activity in nuclear physics (see 
Refs. [7H1I7H1F771F7H1F71I1[R(T1[RT1IR2] ). One one hand, in phe¬ 
nomenological calculations based on relativistic mean field 
models the maximum mass of hyperon stars could still 
reach the 2 M 0 limit (see for instance piHll841IRH] ) on the 
other hand in more realistic calculations based on micro¬ 
scopic nucleon-nucleon interactions, the appearance of hy¬ 
perons is accompanied by a strong softening of the EoS 
which leads to maximum masses much below 2 M 0 and in 
some case even below I. 4 M 0 [5M571IE5j . A possible way 
out is that the hadronic EoS is so stiff that even for the 
2Mq star the central density is below the threshold for 
the appearance of hyperons: a possible example has been 
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Fig. 8. Mass-radius relation for hadronic stars and observa¬ 
tional constraints indicating the existence of very compact stel¬ 
lar objects. 


given in [53] where at a mass 2.O9M0 the central density 
is of 3.5no and hyperons are not yet formed. This scenario 
is realized if the three-body hyperon- nucleon interaction 
is sufficiently repulsive. Another possible way to add re¬ 
pulsion between baryons is the multi-Pomeron exchange 
potential proposed in jSD] which, again, would allow the 
existence of massive hadronic stars. Unfortunately from 
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Fig. 9. Relation between gravitational mass and baryonic mass 
for hadronic stars and quark stars. At fixed baryonic mass the 
difference between the gravitational mass of a hadronic star 
and a quark star is of the order of O.ISM©. In the insert we 
display also the curve corresponding to nucleonic stars within 
the SFHo EoS (green line). The slightly larger binding energy 
obtained when adding deltas and hyperons allow to fulfill the 
constraint of [70] (blue box). 

the theoretical side one cannot draw a firm conclusion on 
the composition of the 2Mq stars. 

The crucial quantity to measure is clearly the radius of 
compact stars. As already discussed, the presently avail¬ 
able analysis are still affected by large systematic errors. 
While there are some hints in favor of the two families sce¬ 
nario (some analysis suggesting the existence of stars with 
radii of the order of 10km and some other analysis which 
infer radii of about 14km) one cannot yet claim that these 
measurements have found compelling evidence for the two 
families scenario. However, the Neutron star Interior Com¬ 
position ExploreR (NICER) instrument, to be launched 
in 2016, will allow to achieve 5% precision in neutron-star 
radius through rotation-resolved spectroscopy (see for in¬ 
stance EH). Still another possibility is based on the fu¬ 
ture X-ray mission ATHENA-I- which, combined with ac¬ 
curate distance measurements provided by the GAIA ex¬ 
periment, will also allow to obtain information on masses 
and radii. Another possible way of extracting the radii is 
by analyzing the gravitational-wave emission during the 
merger of two compact stars. In particular, the frequency 
of the dominant oscillation mode of the post-merger rem¬ 
nant is directly related to the radius of the non-rotating 
compact star m- Again the expected error is of the order 
of a few hundred meters. With such a precision one would 
possibly finally establish what the internal composition of 
compact stars is. Let us consider the canonical I. 4 M 0 ; 
one can imagine three possible outcomes for the measure¬ 
ment of its radius: i) radius larger than about 13.5 km; ii) 
radius between 11.5 and 13.5 km; iii) radius smaller than 
about 11.5km. These numbers are also suggested by the 
model-independent study of Ref. |7| . 


In the first case, the large radius implies that the EoS 
is stiff. Therefore the central density of the 1.4M0 star 
is small and it is possible that even for the most massive 
stars the central density is smaller than the threshold of 
formation of hyperons and deltas. Only “normal” neutron 
stars would exist. In the second case, the interpretation 
would be much more complicated: several equations of 
state have been proposed which can reach the 2Mq limit 
and predict a radius of about 12km for the canonical neu¬ 
tron star mass. Hybrid stars and hyperonic stars would 
be both possible (in particular the scenario i) discussed in 
the introduction would be favored). In the last case, the 
only available possibility is the two families scenario that 
we are here proposing. Indeed there are no examples of 
equations of state in the literature that can at the same 
time fulfill the 2Mq and predict such small radii for the 
canonical neutron star 0. Within a few years one could 
then expect the NICER experiment to give a final answer 
to the question about the existence of the two families. 

The scenario of two families of compact stars is based 
on the existence of quark stars together with hadronic 
stars. In turn the existence of quark stars implies that the 
Bodmer-Witten hypothesis is correct and therefore also 
small nuggets of quark matter, the so called strangelets, 
must exist. We know very little about strangelets: they 
are expected to be positively charged and to have a very 
small charge to mass ratio, see m- Their mass spectrum 
is instead completely unknown. As a matter of fact, to 
date, there is no experimental evidence of the existence 
of strangelets. A possible way to produce strangelets on 
earth is through heavy ions experiments: the search for 
strangelets at RHIC has produced upper limits of few 10“® 
to 10“^ per central Au-|-Au collision for strangelets with 
mass larger than 30 GeV [SSj. Notice however that the 
production of strangelets in heavy ions collision requires a 
net baryon excess [96] which is unlikely to be obtained at 
RHIC energies. In the context of cosmology a very inter¬ 
esting hypothesis is that strangelets could represent a type 
of “macro dark matter” which would have been produced 
by some post-inflationary process m- In this scenario, 
no “beyond standard model physics” would be required 
to explain the existence of baryonic dark matter. In astro¬ 
physics, strangelets could represent an important compo¬ 
nent of cosmic rays. The events that would most probably 
produce strangelets are the merger of compact stars with 
at least one of the two stars being a quark star, see [55] 
for a discussion on the fragmentation of quark matter into 
strangelets. The search for strangelets in the lunar soil has 
provided no evidence of them from A=42 to A=70 and for 
nuclear charges of 5, 6, 8, 9 and 11 [Slj. New limits on the 


^ In Ref. [^, some examples of non-relativistic microscopic 
calculations of nucleonic matter are shown that fulfill the 2Mq 
and, at the same time, predict the existence of very compact 
stars. There are two problems associated with such equations of 
state and in particular with WFFl m- already for the canon¬ 
ical mass the central density is of about 4no. At such densities 
heavy baryons should be included while in WFFl they are sim¬ 
ply neglected. Second: at densities of about lfm“® it violates 
causality. 
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mass and the flux of cosmic strangelets are expected to be 
available in the near future thanks to the AMS-02 exper¬ 
iment. 

From the theoretical side, to date, there has been only 
one detailed simulation of the merger of two quark stars 
[100] while the possibility of neutron star - quark star and 
black hole - quark star merger have not yet been consid¬ 
ered. An unexpected result of these simulations is that 
in many cases, after the merger, a prompt collapse to a 
black hole occurs and basically no quark matter is ejected. 
In particular, this occurs for values of the total mass of 
the merger larger than about 2.5 — SMq. In our scenario 
quark stars have masses larger than about 1.35Mq and it 
is rather difficult to avoid a prompt collapse. It is therefore 
possible that even if quark matter is absolutely stable the 
flux of strangelets is vanishingly small and not all compact 
stars convert into quark stars as it would result if the cos¬ 
mic strangelets pollution would be significant |101] . The 
two families of compact stars could indeed coexist. 

A.D. would like to thank V. Hislop for the moral sup¬ 
port during the preparation of the paper. G.P. acknowl¬ 
edges financial support from the Italian Ministry of Re¬ 
search through the program Rita Levi Montalcini . 
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